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Abstract: The data collected by the Pierre Auger Observatory are analyzed to search for coincidences 
between the arrival directions of high-energy cosmic rays and the positions in the sky of astrophysical 
transients. Special attention is directed towards gamma ray observations recorded by NASA's Swift 
mission, which have an angular resolution similar to that of the Auger surface detectors. In particular, 
we check our data for evidence of a signal associated with the giant flare that came from the soft gamma 
repeater 1806-20 on December 27, 2004. 



Introduction 

It has long been known that two of the highest 
energy cosmic rays came from directions that are 
within the error boxes of two remarkable gamma 
ray bursts (GRBs) detected by BATSE, with a de- 
lay of roughly 10 months after the bursts [1]. How- 
ever, all searches for coincidences between the ar- 
rival direction of ultra-high energy cosmic rays 
(UHECRs) and GRBs from the Third BATSE cat- 
alogue yield a negative result [2]. These searches, 
of course, may have been distorted by the poor 
angular resolution (about 3 degrees) of the GRB 
measurements. A sensitive anisotropy analysis is 
now feasible using data collected by the Swift mis- 
sion [3] and the Pierre Auger Observatory [4]. In 
the first part of this paper we present the results of 
such an analysis. 

In the second part of the paper, we search for a sig- 
nal from the direction of the soft gamma repeater 
(SGR) 1806-20, which on December 27, 2004 
emitted a hyperflare that saturated many satellite 
7-ray detectors [5]. This unique flare lasted about 
5 minutes (the duration of the initial spike was ~ 
0.2 s), had a peak luminosity of ~ 2 x 10 47 erg/s, 
and a total energy emission of ~ 5 x 10 46 erg [6]. 
The model proposed is that of a "magnetar" (i.e., 
a neutron star with a huge magnetic field, B ~ 
10 15 G) located on the far side of our Galaxy (at 



a distance rs 15 kpc [7]). The origin of the flare 
can be explained as global crustal fractures due to 
_B-field rearrangements liberating a high flux of Ar- 
rays and 7-rays [8]. The exceptional energetics 
of this hyperflare makes of SGR 1806-20 an at- 
tractive candidate source of UHECRs, high energy 
neutrinos, and gravitational waves [9]. Given that 
the source is in our Galaxy, if it were to generate 
high energy neutrons a significant fraction of them 
could arrive at Earth before decaying and would 
point back to the source. Moreover, neutrons that 
decay in flight would produce antineutrinos that in- 
herit directionality [10]. Searches for neutrino and 
gravitational wave emission have been reported by 
the AMANDA [11], AURIGA [12], and LIGO Sci- 
entific [13] collaborations. In all these searches the 
data revealed no significant signals. 

Cosmic Ray Data Sample 

In our analysis we use data collected with the 
Pierre Auger Observatory, located in Argentina 
at Southern latitude 35.2° and Western longitude 
69.3°. We consider events with zenith angle 9 < 
#max = 60°, detected from January 1, 2004 to 
April 1, 2007. We employ the reconstruction 
procedure discussed elsewhere [14]. A total of 
609,161 CRs have passed the selection criteria. 
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Figure 1 : The dots indicate the positions of 284 well-localized GRBs (equatorial coordinates J2000) ob- 
served from January 1, 2004, to April 1, 2007. 192 of them are within the f.o.v. of Auger with a zenith 
angle < 60°, and 234 with a zenith angle < 80°. The boundary of Auger f.o.v. for 6* max = 60° is defined 
by declination = 24.8°. The position of SGR 1806-20 is indicated by an open square. 



UHECRs and GRBs 

A catalogue of 284 GRBs observed with an accu- 
racy of better than 1° (from Januray 1, 2004 to 
April 1 , 2007) was compiled using data primarily 
from the Swift mission complemented by measure- 
ments from additional GRB observing satellites, 
including HETE, INTEGRAL, and others. Out of 
the total GRB sample, 192 bursts are within the 
field of view (f.o.v.) of Auger and only 62 were 
in the Auger f.o.v. at the time of their bursts, i.e. 
&GKB < #max- The GRB sky distribution is given 
in Fig. 1. As expected [15], they cover the whole 
sky isotropically. 

We bin CR events coming from directions defined 
by spherical caps of radii ?9 sop = 5° and ^ scp = 
30° around each GRB position. We determine the 
total number of coincidence candidates by count- 
ing the number of CRs found within each of the 
specified cones. The number of coincidence can- 
didates for 30° is AA 30 o = 1, 980, 577 whereas for 
5°, J\f 5 c = 57, 053. Note that for 30° the number 
of coincidence candidates is larger than the total 
number of CRs in the sample, since a given CR 



may lie within 30° of more than one GRB. Differ- 
ences between the observation times of the GRBs 
and the arrival times of CR events in the same an- 
gular bin were determined. In Fig. 2, we show the 
rates of CR events as function of the GRB-CR time 
difference for the two angular radii i9 sop around the 
positions of GRBs. We consider a 100-day period 
before and after the GRB observation. The his- 
togram corresponds to a 30° cut and the points cor- 
responds to a 5° cone. No significant excess after 
the time of the bursts is evident in the data. 

UHECRs from SGR 1806-20? 

The giant flare of December 27, 2004 from SGR 
1806-20 represents one of the most intriguing 
events captured in almost three decades of mon- 
itoring the 7-ray sky. Such an energetic event 
clearly constitutes a potential candidate for accel- 
eration of UHECRs. Secondary neutrons can be 
produced in collisions of relativistic protons (and 
nuclei) with the ambient plasma. Interestingly, 
those produced with an energy E > 10 18 eV have 
a boosted cr„ sufficiently large to serve as Galac- 
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Figure 2: Rates of CR events as a function of the 
difference between the GRB time and the CR ar- 
rival time. Data falling within 30° of a GRB are 
indicated by the histogram and within 5° by the 
points. For clarity, statistical errors are only shown 
for the 5°-distribution. 
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Figure 3: Rates of CR events from the direction 
of SGR 1806-20 as a function of the 7-CR time 
difference. Conventions are the same as in Fig. 2. 



tic messengers. The decay mean free path of a 
neutron is cT n T n = 9.15(£^/10 18 eV) kpc, the 
lifetime being boosted from its rest frame value, 
t„ = 886 s, to its lab value by r„ = E/m n . Be- 
cause of the exponential depletion, about 20% of 
the neutrons survive the trip at 10 18 eV, and about 
58% at 10 18 - 5 eV. 

The location of the source, right ascension 
18h 08m 39.34s and declination -20° 24' 39.7", 
is within the f.o.v. of Auger, and below 60° for 
about 9 hours per day. At the flaring emission 
21:30:26.5 UTC, its zenith angle was #sgr — 
70.3°, and it remained for a 50 minute interval 
above the horizon. Unfortunately, this is outside 
the currently best understood region of the detec- 
tor, i.e. 6> max = 60°. 

We have repeated the analysis described in the pre- 
vious section for this exceptional burst. In Fig. 3 
we show the results of such an analysis, indicating 
that no significant excess in the CR flux is evident 
after the burst. (The number of events coming from 
the direction of SGR 1806-20 within a 30° cone is 
Af 30 o = 5, 596 whereas for 5°, Af s <, = 139.) 

By extending our data analysis to higher zenith an- 
gles, we have verified that no events have been ob- 
served within a 5° cone during the T — 300 s 
of the flare, where #sgr ~ 70°. The abscence 
of a signal can then be exploited to place an up- 
per bound on the primary neutron flux, without 
assumptions on the Galactic magnetic field. To 
do so, we must determine the effective detection 
area A and the trigger efficiency e(E). We adopt 
the elementary hexagonal cell approach discussed 
in Ref. [16]. On December 27, 2004 an aver- 
age of 364 Auger hexagones were fully active, 
each counting for 1.95 km 2 on the ground (which 
amounts to 0.66 km 2 as seen with an angle of 
w 70°). Hence, the experiment presented A ~ 
239 km 2 to the potential "beam" of neutrons em- 
anating from the source. The trigger efficiency for 
70° is shown in Fig. 4. 

Using the diffuse flux of cosmic rays we estimate a 
background 1, thus Poisson statistics implies an 
upper bound of 3.09 events at 95% CL from neu- 
tron fluxes, d$/dE [17]. Equivalently, for some 
interval A, 



T A 



d® 

dE — e(E) < 3.09 . 
A dE 



(1) 



Search for coincidences between Auger data and astrophysical transients 




17 17.25 17.5 17.75 18 18.25 18.5 18.75 19 19.25 19.5 

Log, (E/eV) 



Figure 4: The rising dot-dashed line indicates the 
trigger efficiency of Auger for 9 = 70°. The 
arrowed circles indicate upper limits on the en- 
ergy weighted flux of neutrons from SGR 1806-20 
(valid in a logarithmic interval A = 1). 

In a logarithmic interval A where a single power 
law approximation (for the integrand) is valid, we 
obtain [18] 



Eq 



dE 



< 



3.09 



E 



T A e(E ) 



(2) 



where Eq is the energy at the center of the loga- 
rithmic interval and we have taken A = 1 (corre- 
sponding to one e-folding of energy). The 95%CL 
upper limits on the energy weighted flux of neu- 
trons are shown in Fig. 4. 

Conclusion 

We used the Auger data sample to search for UHE- 
CRs which are consistent with the position and 
time of astrophysical transients. No such coinci- 
dences were found in the data. 
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